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Nature, in the form of dissipation, inevitably intervenes in our efforts to control 
a quantum system. In this talk we show that although we cannot, in general, 
compensate for dissipation by coherent control of the system, such effects are not 
always counterproductive; for example, the transformation from a thermal (mixed) 
state to a cold condensed (pure state) can only be achieved by non-unitary effects 
such as population and phase relaxation. 



1 Representation of quantum states for dissipative systems 

In closed-system, pure-state quantum mechanics the state of the system is 
usually represented by a wavefunction | , which is an element of the Hilbert 
space 7i. For open quantum systems, however, a quantum statistical mechan- 
ics formulation is necessary since dissipative effects due to the interaction of 
the system with its environment convert pure states into statistical ensembles 
and vice versa. The state of the system must therefore be represented by a 
density operator p, i.e., a positive trace-one operator acting on 7i. It is con- 
venient to expand this density operator in terms of a complete orthonormal 
set of energy eigenstates {\n) : 1, 2, . . . , iV = dimTi} of the system: 
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such that the diagonal elements p„„ in this expansion determine the pop- 
ulations of the energy eigenstates \n), while the off-diagonal elements pnm 
for n ^ m determine the coherences between energy eigenstates. The latter 
distinguish coherent superposition states |^) = J2n=i'^n\n) from statistical 
ensembles of energy eigenstates (i.e., mixed states) p = X]n=i ""^"k*) ("■!• To 
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see the difference between the two, consider pi — ^ 

pi is diagonal and represents a two-level system in an incoherent mixed 
state with equal populations in states |1) and |2) but no correlation between 
both states. Note that such a state cannoi be represented by a wavefunc- 
tion. p2 has off-diagonal elements and diagonalization shows that it rep- 
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resents the coherent superposition state j^*) = + |2)] since we have 

/5=l*)(*l = ^[l,l]^x-i5[l,l]. 

2 Dynamics of dissipative quantum control systems 

For Hamiltonian systems the evolution of the state p{t) with p(io) — Po is 
governed by 

m = mpoUit)\ (2) 

where U(t) is the evolution operator satisfying the Schrodinger equation 

ihjUit) = Hmit), U{0) = /, (3) 
and / is the identity. p(t) also satisfies the quantum Liouville equation 

ihj^m = [H{f),m] = H{f)m - pmit). (4) 

H{i) is the total system Hamiltonian, which depends on the control fields fm- 

M 

H{i)^H0+J2frnit)H^, (5) 

where Hq is the internal Hamiltonian and Hm is the interaction Hamiltonian 
for the field /m for 1 < m < M. 

The advantage of the Liouville equation over the unitary evolution 
equation (|2Jl is that is can easily be adapted for dissipative systems by adding 
a dissipation (super)operator L£)[/5(t)]: 

M 

ihp{t) = [Ho,m] + J2 fmit)[H,n,m]+ihLD[m]- (6) 

m— 1 

Under certain assumptions (semi-group dynamics, norm continuity and con- 
servation of probability), it can be shown that the dissipation operator has 
the form^ 

s 

where the Vs are arbitrary N x N matrices, i.e., bounded operators on the 
Hilbert space H. 
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3 Population relaixation and phase decoherence 

In general, uncontrollable interactions of the system with its environment 
lead to two types of dissipation: population relaxation (decay) and phase 
decoherence (or dephasing). 

The former occurs, for instance, when a quantum particle in state |n) 
spontaneously emits a photon and decays to another quantum state \k), which 
changes the populations according to 

Pmi(0 = -^{[H{i),p{t)])nn + ^ [lnkPkk{t) ~ 7fc„p,m(i)] (8) 

where jknPnn is the population loss for level \n) due to transitions \n) — > \k), 
and 'jnkPkk is the population gain caused by transitions \k) \n). The total 
population relaxation rate jkn is determined by the lifetime of the state |n), 
and for multiple decay pathways, the relative probability for the transition 
|n)^|fc). 

The latter occurs when the interaction with the enviroment destroys the 
phase correlations between states, which leads to a decay of the off-diagonal 
elements of the density matrix: 

Pknit) = ~^{[H{i),p{t)])kn - TknPknit) (9) 

where Tj^n (for k ^ n) is the dephasing rate between |fc) and Note 
that population relaxation always induces dephasing since decay destroys the 
phase correlations between states. However, there may be other sources that 
contribute to the loss of coherence of the system. 

The effects of dephasing and population relaxation can be accounted 
for by adding a dissipation super-operator defined by (L£)[/5(t)])£m = 
iJ-'D)im,nkPnk{t) to the LiouviUe equation. The latter can be represented 
by an N"^ x A''^ matrix whose non-zero elements are 

(J^D)kn.kn — —^kn, {^D)nn.kk = +7«A:: k ^ fl ^ ^ 

nn,nn — y^i„-Af~ "fkn ■ 

Population decay and dephasing allow us to overcome kinematical con- 
straints such as unitary evolution to create statistical ensembles from pure 
states, and pure states from statistical ensembles, which is important for 
many applications such as optical pumping. However, there are instances 
when this is not desirable such as in quantum computing, where these effects 
destroy quantum information. Hence, there are situations when we would like 
to prevent decay and dephasing. A cursory glance at the quantum LiouviUe 
equation for coherently driven, dissipative systems ^ suggests that it might 
be possible to prevent population and phase relaxation by applying suitable 
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control fields such that 

M 

E 

m=l 



fm{t)[H^,p{t)] + IKLdW)] = 0. 



(11) 



Unfortunately, however, a more careful analysis reveals that this is not possi- 
ble, in general, as we shall now show explicitly for a two-level system, or qubit 
in quantum computing parlance. 

4 Dynamics of a dissipative, coherently driven 2-level system 

The Hamiltonian for a driven two-level system with energy levels E\ < E2 is 



Him] =Ho + fi{t)Hi + f2{t)H2 



(12) 



where Ho is the internal Hamiltonian and Hi and H2 are Hamiltonians de- 
scribing the interaction with independent (real- valued) control fields /i {t) and 
/2(t). For a spin system, for example, the two control fields might be two or- 
thogonal polarizations of an electromagnetic field affecting rotations about 
two orthogonal axes. In general, we can assume without loss of generality 
that the internal and interaction Hamiltonians have the form: 



Ha = 
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, H2 = d2 
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where di , d2 are the (real- valued) dipole moments for the transition and u) = 
{E2 — Ei)/h is the transition frequency. 

We can re-write the Liouville equation in matrix form in a higher dimen- 
sional space (Liouville space). Straightforward computation shows 



-|p(t))) = -i[Lo + fi{t)U + h{t)U + iLz5]|p(t))) 



(13) 



where \p{t))) = [pii{t), pi2{t), p2i{t), p22{t)Y dissipative evolu- 

tion is given by 



Ln = a; 
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with 712 being the rate of population relaxation from |2) to |1), 721 the rate of 
population relaxation from |1) to |2), and F the dephasing rate. The control 
action is given by 
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Observe that the matrix elements of the control Liouville operators Li and L2 
are zero where the matrix elements of the dissipation operator hjj are non- 
zero, and vice versa. Thus, no matter how we choose the control fields, we 
cannot cancel the efi'ect of the dissipative terms without introducing additional 
terms such as measurements and feedback 1^ or the ability to control the 
coupling of the system to the reservoir'"'. 

Note that the dissipation operator as defined in (|14|l is equivalent to the 
Lindblad form since inserting 

with r = r — ^(712 + 721) into (TJl gives 

_ — 721PII —^721/012 
-5721/021 721PII 

^ /^-72ipii + 712P22 -rpi2 
V -rp2i 721P11 - 712P22 

which is equivalent to Ijo\p{t))) with ho as in |(7J). 
5 Dissipation and entropy 

One of the main consequences of dissipation is that interactions of the system 
with a bath (environment) can change the entropy and purity of the system. A 
useful measure of the entropy and purity for our purposes is 1— Tr[/5(i)^], which 
essentially determines the Renyi entropy of the system, although the latter 
is often defined to be — logTr[p(i)^]. For a non-dissipative, coherently driven 
quantum system the entropy is conserved because the evolution must remain 
unitary. Dissipation provides new opportunities for control by enabling us to 
reach states outside the orbit of the initial state under the unitary group El 
especially states whose entropy differs from the initial state. Dephasing, for 
instance, converts coherent superposition states into uncorrelated statistical 
mixtures of energy eigenstates and hence enables us, in principle, to convert 
a given pure state into an arbitrary mixed state by creating a superposition 
state using coherent control, and letting the coherences decay. Population 
relaxation allows us to convert a high entropy mixed state into a (zero entropy) 
pure state and vice versa. 

5. 1 Conversion of a pure state into a mixed state 

When only pure dephasing occurs a coherent superposition of energy eigen- 
states l^f) = J2n=i'^n\'ri) with X]n=i CnC* — 1 dccays into a statistical mix- 
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ture of the states \n) with a discrete probability distribution Wn = |c„p for 
1 < n < A'^. For dephasing times much greater than the control time we can 
design a control field that transforms the initial state into the required coher- 
ent superposition without worrying about dephasing, and then turn the field 
off to let dephasing transform this superposition state into the desired mixed 
state. However, if significant dephasing occurs during the coherent control 
phase, either due to rapid dephasing, or because the control process takes too 
long, then this approach will fail. 

For instance, consider a system with two non-degenerate energy levels. 
Suppose we wish to transform the pure state |1) into a statistical mixture of 
the states |1) and |2). Based on geometric control for non-dissipative systems, 
we might try to apply a resonant Gaussian control pulse with effective pulse 
area ^, which would create the superposition state j^") — + |2)) in 

the non-dissipative case, and hope that decoherence will convert this state 
into the desired mixed state. Our control calculations indicate, however, that 
this scheme will fail for dephasing rates of the order of the Rabi frequency of 
the control pulse. However, straight-forward optimization with respect to the 
effective pulse area and length of the control pulse indicates that the pulse 
length and pulse area can be chosen such as to achieve the desired result. For 
instance, by increasing the effective pulse area of a Gaussian pulse lasting 50 
vibrational periods from the predicted value of ^ for a closed system to O.SItt, 
we were able to create the desired maximum entropy state for a dephasing 
rate F = 0.1 in just over 50 vibrational periods. 

5.2 Conversion of a mixed state into a pure state 

A perhaps even more important application of controlled dissipative dynamics 
in quantum optics is optical pumping to drive a mixed-state system into a 
desired pure state using a combination of coherent control and population 
relaxation from an excited state. For instance, suppose we have a cloud of cold 
atoms whose electronic ground state is three-fold degenerate. If the system 
is not prepared in a particular pure state, it will usually be in a mixture of 
the three degenerate substates, which we shall denote by |2) and |3) for 
simplicity. For many applications, e.g., in quantum computing, it is crucial to 
prepare the system in a certain pure initial state. As we have seen, this is an 
aim generally impossible to realize by coherent control alone. To be able to 
take advantage of spontaneous emission to increase the purity of the system, 
we must couple the ground state to an excited electronic state with a finite 
lifetime. The sublevels of the ground and excited states can be coupled in 
various ways depending on the polarization of the field. The trick is to select 
the right coupling. 

For example, suppose the upper level is also three-fold degenerate and 
the coupling induced by the control field is as indicated in figure ^ i.e., the 
field couples states |2) and |5), as well as |3) and |6). The excited states can 
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emit a photon and return to one of the ground states. Certain transitions are 
prohibited by atomic selection rules; the allowed decay modes are indicated 
in figure n (right). The simplest optical pumping schemes involve applying 
a constant amplitude field resonant with the transition frequency between 
the two levels and suitably polarized to couple only the levels indicated in 
figure n Without population relaxation due to spontaneous emission, the 
field merely leads to population cycling between states |2), |5), and |3), |6), 
respectively. Adding population relaxation changes the effect of the control 
field dramatically, leading to an accumulation of the population in state |1) 
as figure 13 shows. If the control field is applied for a sufficiently long time, all 
the population will eventually accumulate in state 

In the previous optical pumping scheme a simple constant amplitude res- 
onant control field was sufficient to achieve the objective of driving the system 
into the desired pure state. However, this is not always the case. Some appli- 
cations of optical pumping such as laser cooling of internal molecular degrees 
of freedom rely on the interplay of carefully selected control pulses and dis- 
sipation. For example, a molecular vapor at room temperature consists of a 
statistical mixture of molecules in many different ro- vibrational states. Due to 
many closely spaced energy levels and lack of selection rules, there are many 
possible transitions with various transition probabilities that can be excited 
by applying a control field, and many different decay pathways. The situation 
is further complicated by the fact that the timescales for coherent control and 
population relaxation are often quite different. The problem thus appears to 
be nearly hopeless. Yet, it has been shown that this problem can be addressed 
successfully using optimal control for dissipative systems and creative control 
strategies 1^. 

An approach that is especially promising for systems where the timescales 
for control and dissipation are quite different (as in our molecular cooling 
problem) involves breaking the problem up into a sequence of excitation and 
relaxation steps. The goal in each step is to use control theory to design 
control fields to transfer the system from its initial state to a kinematically 
equivalent, dynamically reachable state, which has the same entropy but is 
likely to decay into a state with lower entropy. In principle, the entropy of 
the system can be decreased until it is zero and the system is the desired pure 
state. The main difficulty of this approach is the choice of suitable target 
states for each optimization step, which requires a good understanding of 
the effects of population relaxation and dephasing on various kinematically 
equivalent states, in order to assure that the selected states will decay into a 
lower entropy state. 

6 Conclusion 

We have shown that although coherent control of the system dynamics in 
absence of feedback or the ability to alter the interaction of the system with its 
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Figure 1. Optical pumping for a degenerate two- level system. Transition diagrams for the 
control field and population relaxation 
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Figure 2. Optical pumping for a degenerate two-level system. Without population relax- 
ation, the coherent control field induces population oscillations between states |2) and |5} 
as well as |3) and |6) (Rabi oscillations) and the entropy of the system remains constant 
(left). Population relaxation dramatically changes the effect of the control field, leading to 
an effective pumping of population into the lower left sublevel and an entropy reduction 
as the system approaches a pure state (right). 



environment cannot compensate for dissipative effects sucli as depfiasing and 
population relaxation in open systems, such effects need not be detrimental 
to control, but are in fact crucial for many interesting applications such as 
ensemble preparation or system purification by laser cooling. 
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